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53 ■ ABSTRACT 



Based on the gamma-ray burst (GRB) event rate at redshifts of 4 < z < 
12, which is assessed by the spectral peak energy-to-luminosity relation recently 
found by Yonetoku et al., we observationally derive the star formation rate (SFR) 
for Pop III stars in a high redshift universe. As a result, we find that Pop III 
stars could form continuously at 4 < z < 12. Using the derived Pop III SFR, we 
attempt to estimate the ultraviolet (UV) photon emission rate at 7 < z < 12 in 
which redshift range no observational information has been hitherto obtained on 
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ionizing radiation intensity. We find that the UV emissivity at 7 < z < 12 can 
make a noticeable contribution to the early reionization. The maximal emissivity 
is higher than the level required to keep ionizing the intergalactic matter at 
7 < z < 12. However, if the escape fraction of ionizing photons from Pop III 
objects is smaller than 10%, then the IGM can be neutralized at some redshift, 
which may lead to the double reionization. As for the enrichment, the ejection of 
all metals synthesized in Pop III objects is marginally consistent with the IGM 
metallicity, although the confinement of metals in Pop III objects can reduce the 
enrichment significantly. 

Subject headings: early universe — gamma-ray: bursts — stars: formation 

1. Introduction 

The recent observation on the cosmic microwave background by Wilkinson Microwave 
Anisotropy Probe ( WMAP) has revealed that the dark age in the universe ended quite early 
by the reionization at z r ~ 20+g (Bennett et al. 2003; Spergel et al. 2003). But, the origin 
of ionizing radiation in such an early epoch and the subsequent ionization history after the 
early reionization are not elucidated very well. The first possibility for the ionizing source is 
Pop III stars (Cen 2003a; Ciardi Ferrara, & White 2003; Wyithe & Loeb 2003; Somerville & 
Livio 2003; Fukugita & Kawasaki 2003), and the other possibilities are mini quasars (Madau 
et al. 2004) or BH accretion (Sasaki & Umemura 1996; Ricotti & Ostriker 2004b). The early 
reionization by Pop III stars is recently under debate. Sokasian et al. (2004) argue that 
the reionization by Pop III stars is marginally compatible with the Thomson optical depth 
measured by the WMAP, r e = 0.17±0.04, even if the escape fraction of ionizing photons from 
early formed objects is unity. In addition, Ricotti & Ostriker (2004a) explore the feedback 
effect by supernova (SN) energy input in Pop III objects, and find that the suppression of Pop 
III star formation by SNe is so intensive that Pop III stars cannot contribute significantly to 
the cosmic reionization. The strong suppression of Pop III star formation or the early shift 
to Pop II or Pop I star formation may bring the vanishment of ionizing radiation, resulting in 
the re-neutralization of intergalactic medium (IGM). This may lead to the so-called double 
reionization (Cen 2003b; Wyithe & Loeb 2004). 

Also, Pop III SNe can be the origin of IGM heavy elements (Nakamura & Umemura 
2001; Scannapieco, Ferrara, & Madau 2002; Ricotti & Ostriker 2004a). The detailed analysis 
on quasar metal absorption lines have shown that the IGM is enriched by heavy elements 
from a level of ~ 10~ 3 Z Q at z < 4.3 to ~ 10~ 4 Z Q at z ~ 5.3, where Z Q is the Solar 
abundance (Cowie and Songaila 1998; Songaila 2001). It is quite perplexing that heavy 
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elements are observed commonly and uniformly in the low density space where no bright 
galaxies are found. Although the outflow from small proto-galaxies could be responsible 
for the enrichment of IGM (Madau, Ferrara, & Rees 2001; Mori, Ferrara, & Madau 2002), 
the much earlier enrichment by Pop III objects may provide a potential solution for the 
uniformly distributed metagalactic heavy elements. Hence, Pop III star formation history in 
high redshifts could be a key for the early enrichment. 

We cannot observe directly Pop III stars at 7 < z < z r , but GRBs at high redshifts can 
be a potential tool to probe the formation history of Pop III stars. Since the discovery of 
X-ray afterglow from GRB 970228 (Costa et al. 1997; van Paradijs et al. 1997), the distances 
to bright GRBs are known to be cosmological (Metzger et al. 1997; Greiner 2004). The 
most distant GRB hitherto identified spectroscopically is GRB 000131 at z — 4.5 (Andersen 
et al. 2000; Greiner 2004). However, there are numerous much weaker GRBs without known 
distance, listed in the archive by the Burst And Transient Source Experiment (BATSE) on- 
board the Compton gamma-ray astronomy satellite (BATSE 2004). Recently, Yonetoku 
et al. (2004) estimated distances of 689 BATSE weak GRBs, based on the spectral peak 
energy-to-luminosity relation (Amati et al. 2002; Yonetoku et al. 2004), and found that most 
of GRBs are more remote up to z ~ 12. A recent great advance on GRBs is the finding of 
the obvious association between GRB 030329 and SN 2003dh (Hjorth et al. 2003; Greiner et 
al. 2003; Matheson et al. 2003; Price et al. 2003; Uemura, M. et al. 2003). This association 
firmly established the connection between a GRB and an energetic supernova (hypernova), 
which is the collapse of a massive star (Woosley et al. 1999; Paczynski 1998). Since Pop 
III stars are likely to form in a top-heavy initial mass function (IMF) (Bromm, Coppi & 
Larson 1999; Abel, Bryan, & Norman 2000, 2002; Nakamura & Umemura 2001), a portion 
of GRBs at high redshifts possibly result from Pop III hypernovae. Actually, Heger et al. 
(2003) estimated the fraction of Pop III GRBs using the IMF by Nakamura & Umemura, 
and found that several percent of Pop III stars result in GRBs. 

In this Letter, we derive the star formation rate (SFR) of Pop III stars observationally 
for the first time out to z ~ 12 by use of the absolute GRB event rate. Then, using the Pop 
III SFR, we assess the UV photon emission rate at 7 < z < 12, in which redshift range no 
observational constraint has been hitherto placed on ionizing radiation intensity. Also, we 
analyze the early metal enrichment of IGM by the derived Pop III SFR. Here, we assume a 
ACDM universe with fi A = 0.7, Vt m = 0.3 and H = 72 kms^Mpc- 1 (Bennett et al. 2003; 
Spergel et al. 2003). 
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2. GRB Rate and Pop III SFR 

2.1. Absolute GRB Event Rate 

Yonetoku et al. (2004) published the relative value of the GRB event rate up to z ~ 12, 
based on the spectral peak energy-to-peak luminosity relation. In this Letter, we need to 
estimate an absolute value of GRB event rate (pgrb) in the early universe. In order to do 
that, it is necessary to make corrections for jet-collimation, because the GRB emission is 
collimated in the order of 0.1 radian (Kulkarni et al. 1999; Frail et al. 2001). If the reported 
evolution in luminosity (Yonetoku et al. 2004) can be attributed to the evolution of jet- 
collimation, we can correct the collimation using the observed evolution in luminosity with 
the index k defined as (1+z) 2 ' 6k , where 2.6 represents the observed evolution in luminosity 
(Yonetoku et al. 2004). After the correction against jet-collimation, the absolute event 
rate pgrb in a comoving volume is obtained. It is not quite certain whether the observed 
evolution in luminosity is fully due to the collimation or partially due to intrinsic evolution 
in luminosity. So, we use the intermediate case of /t=0.5 for simplicity in this paper. The 
resultant rate in units of Gpc -3 yr" 1 is shown in Fig. 1 together with the observed data 
with the systematic and statistical errors. The typical rate at z ~ 12 is roughly 10 4 Gpc~ 3 
y r -i f or ^ e case f k=0.5. 

2.2. Pop III SFR 

The absolute GRB rate allows us to estimate the SFR of Pop III stars with the assistance 
of the Pop III IMF. Because of the lack of heavy elements, the IMF for Pop III stars is likely 
to be top heavy with massive stars of > 5OM , in contrast to the present-day stars (Abel, 
Bryan, & Norman 2002; Nakamura & Umemura 2001; Tan & McKee 2004). The possible 
formation of even lower mass Pop III stars with 1 to 50 M is also pointed out (Nakamura 
& Umemura 2001). But these stars might be a quite small contribution to the UV emission 
rate or GRB event rate, and therefore omitted in the present analysis. 

The theory on the stellar evolution of extremely metal poor (< 1O~ 4 Z ) stars (Heger 
et al. 2003) shows that the stars with about 8 to 40 M Q will evolve and explode as type-II 
SNe, leaving a neutron star or a black hole and supplying heavy elements into space. The 
stars with 40 to 140 M will produce UV photons during their evolutionary pass, but will 
not supply any heavy element into space, because they will collapse completely, forming a 
black hole. The stars with 140 to 260 M will explode as a pair-instability SNe, forming a 
black-hole and supplying heavy elements into space. Among these stars, the key interest is 
the stars with 100 to 140 M , which can result in GRBs as progenitors. This mass range 
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is sometimes called an energetic SN (hypernova). Hence, in metal poor environments, GRB 
events occur only in the narrow mass range in the Pop III IMF. 

Here, we employ the Pop III IMF by Nakamura & Umemura (2001), which is in the 
form of a double-peaked function with sharp cutoffs with power-law shape: a — (5 — 1.35 
and cutoffs: m p i = 1M©, and m p2 = 5OM respectively (see Nakamura & Umemura 2001 for 
details). As stated above, we consider only a high- mass part of the IMF. The normalization 
factor can be determined by the GRB event rate, Pgrbj which corresponds to a mass range 
of 100 — 14OM . Then, if the Pop III fraction /q^ b in the observed GRB rate is given, we 
can derive the Pop III SFR in terms of 

p™ ~ 3 x 10 2 /™ B p GRB M Gpc- 3 yr- 1 (1) 

after an integration over the IMF. Actually, /grb i s a function of time. Heger et al. (2003) 
show that the GRB range quickly expands for Z > 10~ A Z & (Pop I /ll), but simultaneously the 
IMF is likely to become the present-day type (Salpeter) IMF for the enriched gas (Bromm, 
Ferrara, Coppi, & Larson 2001). As a result, the mass fraction of GRBs originating from 
Pop I/II stars is roughly twice as high as that of Pop III GRBs. Hence, if J™ is the fraction 
of Pop III stars, the /grb * s gi ven by 

/^b^ 0.5/^/(1-0.5/™). (2) 

It is expected that /J 11 is near unity in an early epoch, and then decreases with the metal 
enrichment. Scannapieco, Schneider, & Ferrara (2003) have derived the /J 11 as a function 
of redshift. Theyhave analyzed the dependence of /J 11 upon the total energy input into 
outflows from Pop III ejecta per unit gas mass, Eg 11 in units of 10 51 ergs M Q -\ The model 
with E l g l = 10~ 4 ' 5 is roughly corresponding to /J 11 = 40 — 50%, while the model with 
£ ni = 10- 3 - 5 is to fl u = 10%. Hence, combining (2) with (1), it is concluded that the 
observed GRB event rate implies the continuous formation of Pop III stars at 4 < z < 12. 



3. Reionization and Enrichment 

3.1. UV Emissivity and Reionization History 

The massive stars with more than 50 M radiate photons in the black body of ~ 10 5 K, 
which are energetic enough to ionize neutral hydrogen. Their luminosities should be near the 
Eddington luminosity (Schaerer 2002), so ~ 10 62 (M/M ) UV photons are emitted during 
their lifetime of ~ 10 6 yrs. For the Pop III SFR assessed above, we can evaluate the total 
UV photon emission rate per unit cosmological comoving volume as 

AT 7 * 3 x 10 62 pi n Gpc-V- 1 (3) 
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for stars heavier than 50 M , where p 1 } 1 is given in units of M Gpc _3 yr _1 . The resultant UV 
emissivity by Pop III stars is shown in Fig. 2. On the other hand, the UV photon emission 
rate required to keep ionizing hydrogen in the universe is estimated as Cnn{z = 0)/t rec 
(Madau, Haardt & Rees 1999), where C is the clumping factor of IGM and t rec is the 
recombination time, t rec = 1100(1 + z)~ 3 (Qbh 2 /0.0224)~ 1 Gyr. Several authors have argued 
the effect of dumpiness on the reionization (Gnedin & Ostriker 1997; Nakamoto et al. 2001; 
Springel & Hernquist 2003). Since the dumpiness shortens the recombination time in HII 
regions, then the clumping factor should be estimated by C = (^Hii)/( n Hii) 2 - When all the 
matter is totally ionized, C becomes equal to the clumping factor of baryons, C = (nl)/(nb) 2 . 
But, if the IGM is partially ionized where density peaks are self-shielded, C is smaller than 
( n b)/( n b) 2 - Therefore, C = 1 corresponds to the minimum level requisite for the IGM 
ionization, while C = (n 2 ) / (rib) 2 does to the maximum level. In Fig. 2, the UV emission 
rate for reionization is shown for C = 1 and C = (n 2 ) / {rib) 2 which is taken from Gnedin & 
Ostriker (1997). The most recent results of WMAP suggest much higher clumping factor than 
their value but we should wait for the second WMAP results for more accurate disccusion. 
As seen in Fig. 2, the UV photon emission rate inferred from the GRB events exceeds the 
maximum level for the IGM ionization at all epochs of 7 < z < 12, for the E™ = 10~ 4 ' 5 
model or for the Pop III GRB fraction larger than several 10%. However, if the escape 
fraction of ionizing photons is much smaller than 10%, the photon emission rate can fall 
below the maximum rate. Especially, when larger systems form at lower redshifts in a CDM 
universe, the escape fraction may be reduced significantly (Kitayama et al. 2004; Whalen 
et al. 2004). Then, the IGM may be neutralized at some redshift, At lower redshifts, Pop 
I/II stars would make more contributions to ionizing photon budget, which may reionize the 
IGM again. 



3.2. Early Enrichment of IGM 

Next, we consider the enrichment by heavy elements synthesized in Pop III objects. 
Since the amount of mass ejection from energetic SNe (hypernovae) and from pair instability 
SNe into space is very roughly ~ 0.1 and ~ 0.5 of the initial mass for 100 — 140 M and 
140 — 26OM respectively (Heger et al. 2003), we can estimate the heavy element production 
rate into the IGM as 0.3p£ n M Gpc~ 3 yr" 1 . Using the Pop III SFR, the metallicity of 
the IGM enriched by Pop III stars is estimated to be 2.2 x 1O~ 4 Z for E™ = 10~ 3 - 5 and 
1.2 x 10- 3 Z Q for Ef 1 = 10~ 4 - 5 . The model with E™ = 10~ 4 - 5 appears to overproduce the 
IGM metals, if compared to ~ 10~ 4 Z Q at z ~ 5.3 (Songaila 2001). However, it should be 
noted that all synthesized metals cannot be spread out from Pop III objects (Norman et al. 
2004). Then, even the model with E^ u = 10~ 4,5 may not provide sufficient metals for the 
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IGM metallicity. 

4. Conclusions and Discussion 

In this paper, we have estimated the absolute GRB event rate out to z ~12 observation- 
ally with the correction for collimation. There have been also several efforts to estimate a 
GRB event rate in the early universe, using the independent method (Murakami et al. 2003; 
Fenimore & Ramirez- Ruiz 2000; Schaefer 2001). All the results are broadly consistent with 
each other and show that there is no break in the GRB event rate as a function of redshift. 

Based on the GRB event rate, we have found that Pop III stars could form continuously 
at 4 < z < 12. This implies that Pop III star formation is not likely to be strongly suppressed 
in early epochs. This is qualitatively consistent with the recent study on the star formation 
history in dwarf spheroidals in the Local group, which shows the continuous star formation 
even after the reionization (Grebel & Gallagher 2004). Also in recent theoretical studies, it is 
shown that dwarf galaxies can form owing to the self-shielding even during the reionization 
(Susa & Umemura 2004), and the star formation can be induced by SN-driven shock (Mori, 
Umemura, & Ferrara 2004). 

As for the reionization, it has been found that the UV emission from Pop III stars can 
make a noticeable contribution to the early reionization. If the escape fraction of ionizing 
photons is large enough, then the UV emission from Pop III stars can keep ionizing the 
IGM at 7 < z < 12. But, if it is smaller, it can lead to the re-neutralization of IGM at 
some redshift, which may result in the double reionization (Cen 2003b; Wyithe & Loeb 
2004). Regarding the enrichment, a high Pop III fraction in high-z GRBs may lead to the 
overproduction of IGM metals if all synthesized metals are spread out into the intergalactic 
space. But, if they are confined around Pop III objects, it may be difficult to account for 
IGM metals solely by Pop III stars. 

Finally, it is noted that there still remain a few less robust factors in this type of 
analysis. The most unknown factor is the validity of the spectral peak energy-to-luminosity 
relation for GRBs out to z ~ 12. Although Yonetoku et al. (2004) derived a firm relation 
out to about z ~ 5 based on the observed redshifts, this relation to higher redshifts is not 
confirmed. The second key is the correction for collimation of GRBs. GRBs are really 
collimated, but we do not know a degree of collimation in luminosities for GRBs at z > 5. 
These should be confirmed by detecting GRBs out to z ~ 10 in the mission like the SWIFT 
GRB satellite. Also, we should check the consistency of the obtained Pop III SFR with the 
infrared background. 
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Fig. 1. — Absolute GRB event rate, Pgrb> per unit cosmological comoving volume, after the 
jet-collimation evolution correction. The emission of GRBs is collimated with typical angle 
of 0.1 radian (Kulkarni et al. 1999; Frail et al. 2001) and probably shows a cosmological 
evolution. We estimate the correction of jet-collimation with a parameter k as (1+z) 2 ' 6 *, 
where 2.6 represents the observed luminosity evolution. The squares with error are the 
observed GRB event rate after the jet-collimation evolution correction with k=0.5. The 
raw-data is used in the paper by Yonetoku et al. (Yonetoku et al. 2004). The dotted lines 
indicate the systematic error of the formation estimate. The three solid lines show the 
differences of the cosmological jet-collimation angle evolution correction with k = 0, 0.5 and 
1.0. The case of k = corresponds to no jet-collimation evolution included. The rough value 
of Pgrb is 10 4 Gpc~ 3 yr _1 at z ~ 12. 
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Fig. 2. — The UV photon emission rate by Pop III stars in a comoving volume inferred from 
the Pop III GRB event rate. A straight line is the emission rate assuming that all GRBs 
originate from Pop III stars, where the raw data are plotted with error bars. Two dashed 
lines show the emission rate using the fractions of Pop III stars based on the models of 
Scannapieco, Schneider, & Ferrara (2003) with E™ = 10 -3 ' 5 and Eg = 10~ 4 ' 5 , where Eg 
is the total energy input into outfows from Pop III ejecta per unit gas mass. The boundary 
of the gray region shows the minimum level required to keep ionizing IGM, assuming the 
clumping factor of C — 1. A solid curve is the maximum level requisite for the IGM 
ionization, assuming C = (nl)/(rib) 2 based on the simulation by Gnedin & Ostriker (1997). 



